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Abstract 

In a previous paper it was shown that the distribution of measured 
values for a retrodictively optimal simultaneous measurement of po- 
sition and momentum is always given by the initial state Husimi 
function. This result is now generalised to retrodictively optimal si- 
multaneous measurements of an arbitrary pair of rotated quadratures 
xg-^ and pg^. It is shown, that given any such measurement, it is pos- 
sible to find another such measurement, informationally equivalent to 
the first, for which the axes defined by the two quadratures are per- 
pendicular. It is further shown that the distribution of measured val- 
ues for such a meaurement belongs to the class of generalised Husimi 
functions most recently discussed by Wiinsche and Buzek. The class 
consists of the subset of Wodkiewicz's operational probability dis- 
tributions for which the filter reference state is a squeezed vaccuum 
state. 
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1. Introduction 

During the last few years there has been considerable progress in the problem 
of simultaneously measuring both the position and the momentum of a quantum 
mechanical system]^, ||. In several recent publications ^, ||] we have discussed 
how to characterise the accuracy of, and disturbance caused by such measurements. 
One approach to the problem is that based on the concept of an "unsharp observ- 
able" ^, 0, H, D- This approach has recently been criticised by Uffink |]lOj. In 
the papers just mentioned we took a rather different approach, based on methods 
developed by Braginsky and Khalili . It appears to us that these methods have 
certain advantages, both conceptually (they clarify what is meant by the term "ac- 
curacy" in a quantum mechanical context ||^), and practically (they facilitate the 
calculations Q). 

Another advantage of these methods is that they give additional insight into 
the physical significance of the Husimi function The fact that the Husimi 

function describes the distribution of measured values for many particular cases of 
joint measurement processes is, of course, well known ^. In ref. (also see 
Prugovecki and Ali 1^) we showed that the Husimi function actually has a much 
stronger, universal property: namely, it gives the distribution of results for any 
retrodictively optimal measurement process {i.e. any process which is retrodictively 
unbiased, and which maximises the retrodictive accuracy). 

The fact that the Husimi function gives the distribution of results whenever the 
measurement is retrodictively optimal, and is otherwise independent of the details 
of the particular process employed, could be interpreted to mean that the Husimi 
function plays the same role for joint measurements of x and p that is played by 
the function | (x | V') P for single measurements of x only. 

The purpose of this paper is to show that a similar universal property holds for 
generalised Husimi functions of the form 

Qgenix,p) j dx'dp' exp [- {a{x - x'f + 2c{x - x'){p~p') + b{p-p'f)] 

X W{x',p') (1) 

where a, b and c are real, a, 5 > 0, a6 — = 1, and where W is the Wigner function. 
These functions are operational distributions of the type defined by Wodkiewicz 
H, |l^ . They are the distributions which result when the filter reference state 
(or "quantum ruler" ) is an arbitrary squeezed vacuum state . They have been 



discussed by Halliwell Wiinsche [IS) and Wiinsche and Buzek [|19[ . 



In ref. 1^ we considered retrodictively optimal measurements of x and p. How- 
ever, this is clearly not the only way to determine the location of a system in phase 
space. What would happen if, instead of determining x and p, one were to make 
a retrodictively optimal measurement of an arbitrary pair of rotated quadratures, 
not necessarily at 90° (see Fig. |l|)? Such measurements are possible (using a suit- 
ably modified form of homodyne detection ||], for example). We will show that the 
outcome of such a measurement is always described by a generalised Husimi func- 
tion of the kind defined by Eq. (0). As with our previous result this is a universal 
statement: it only depends on the measurement being retrodictively optimal, and 
is otherwise independent of the details of the particular process employed. 

The main difficulty in proving this result comes from the fact that the measure- 
ments we consider are characterised by three independent parameters [namely,two 
angles 9,(p to specify the oblique coordinate system (see Fig. |l|), and a parameter 
A to specify the relative accuracy of the measurements of the two quadratures]. 
On the other hand it only needs two parameters to specify a distribution of the 
type defined by Eq. (|l|). It follows, that corresponding to any given distribution. 
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there is an infinite set of informationally equivalent measurements. It is the prob- 
lem of characterising these sets, and giving a precise definition of "informational 
equivalence," which will mainly concern us in the following. 

One might also ask what is the significance of distributions which are like the 
ones considered in this paper in that they are obtained from the Wigner function by 
smoothing it with a Gaussian convolution, but in which the determinant ab — c^ > 1 
(corresponding to an impure filter reference state) ^ [l^, ^ . In ref . Q 

we showed that, in the special case of the Arthurs-Kelly process, such functions 
describe the outcome of measurements producing the smallest possible amount of 
disturbance for a given, sub-optimal degree of accuracy. It is natural to wonder 
whether this property is also universal, and whether it also generalises to the case of 
simultaneous measurements of an arbitrary pair of rotated quadratures. However, 
that is a question which we leave to the future. 

2. Linear Canonical Transformations 

Squeezed states arise as a result of making linear canonical transformations of 
the creation and annihilation operators We begin by describing the parame- 
terisation of these transformations which will be employed in the sequel. 

Consider a system, having one degree of freedom, with position x and momentum 
p. In some applications x, p are dimensionless to begin with. If not they can be 
made dimensionless, by making the replacements x jx, p — > j^p, where I is a 
fixed, in general arbitrary constant having the dimensions of length. In the sequel 
we will always assume that this has been done, so that [x,p] = i. 

We are interested in transformations of the form 

where M belongs to the group SL(2,R) [which is isomorphic to Sp(2,M)]. The 
group may be parameterised as follows. Given any matrix M G SL(2, R) there exist 
unique 9 in the range —'k<9<'k and unique in the range —j<4><j such 
that 

_ / a cos(6' + (j)) a sm{9 + 
\-f3 sm{9 - (/)) p cos{9 -(j))) 

for suitable positive constants a, (3 (see Fig. |^). The requirement that det A/ — 1 
means that a/3 cos 20 — 1. We may therefore write 

for unique A in the range < A < oo. We will refer to 9 as the rotation, as the 
obliquity and A as the resolution. 

The matrix M defines a metric on phase space, with metric tensor M"^ M: 

dslj = dxlj + dpli = (dx dp) M'^M (^^^ (3) 

(M"^ being the transpose of M). For a given system some choices of the matrix M, 
and therefore some choices of metric, will be more natural than others. However, 
if one wants to keep the discussion completely general, so that the detailed nature 
of the system is left unspecified, then one must regard the different choices for M 
as all being on the same footing — corresponding to the well-known fact, that there 
is (in general) no natural metric on phase space. 

From this point of view, the choice of one particular conjugate pair x and p as 
basic, xm and pM being defined in terms of them, must be regarded as arbitrary. 
Assignments of angle also depend on the choice of metric tensor, and must likewise 
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be regarded as arbitrary. It follows, that in a general context (though perhaps not 
in applications to a particular type of system) , no fundamental significance attaches 
to the distinction between oblique axes {(f) ^ 0) and perpendicular axes ((/) = 0). 

It will be shown below that to each M g SL(2,R) there corresponds a retrodic- 
tively optimal measurement; and that two such matrices define the same metric if 
and only if the corresponding measurements are informationally equivalent. 

M has the decomposition 

M^TxS^Re (4) 

where 

/ cos6' sin6'\ „ , Tr-fcoscf) sin(/i\ fj 0\ , . 

Re = { -a a\ ^<t> = V sec 20 . \] = a (5) 

smW cosO J ^ ^ \^sm0 cosqyj \0 \ J 

Define a = -^{x + ip). Let tjg be the unitary rotation operator, and V^, Wx the 
unitary squeeze operators jl^ defined by 

Ue — exp [—i9 a^a] 

% = exp[itanh"^ (tan(/)) {a? + a) '^)] 
M^A =exp[ilnA(a2-at2^j 



Then 



(6) 

n (;) % - g) 

Wl (;) = (I) (8) 



Hence 



3. Informationally Equivalent Measurements 

The purpose of this section is to show how the set of retrodictively optimal 
measurements divides into subsets of informationally equivalent measurements. 

Suppose that we make a retrodictively optimal measurement of the conjugate 
observables xm, Pm, of the kind described in ref. Let /ixM and /ipM be the 
pointer observables giving the results of the measurements of xm and pm respec- 
tively. Let C/mcas be the unitary evolution operator describing the measurement 
interaction, and let exMi, epA/i be the retrodictive error operators 

ixMi = C^mcas I^^mU^^^^ - Xm 
epA/i = t^moas APA/C/mcas ~ -PAf 

as defined in ref. ||^. If the measurement is retrodictively optimal there exists |^ 
fixed T such that 

- T 
((-0 «> (/-apl exA/i F «> (/-ap)) ' = Acia;Af = 

1 I 

((V" ® -^apl epAfi |V' ^ (/lap)) ^ = AciPM 



V2 



T 



for every normalised initial system state lip) (where |0ap) is the initial apparatus 
state, and AciXi\i , AdPM are the maximal rms errors of retrodiction p, M ) . 
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There is no loss of generahty in confining ourselves to balanced measurements, 
for which r = 1. In fact, suppose that r 7^ 1. Define 

It can be seen that measuring the observables xm, Pm to retrodictive accuracies 
±^ and respectively is equivalent to measuring the observables xn, Pn 

both to the same retrodictive accuracy i"^- 

Let M be the matrix with rotation 6, obliquity (j) and resolution A, as in Eq. (||). 
Define 

xe+if, — cos{9 + (j))x + sin(6' + 0)p = — sin{6 — 4>)x + cos(6' — 

Then 



\/scc 2(f> 



a;M = ^ Pa/ = AVsec20pe-0 

It follows, that making a retrodictively optimal, balanced measurement of xm, 
Pm is equivalent to making a retrodictively optimal measurement of a;e+0, Pe-^ 
to accuracies ±-^=^= and ±^^^^^===^ respectively. This equivalence means that 
there is associated, to each retrodictively optimal measurement of a pair xg+^, 
Pe-cf,, a unique matrix £ SL(2,M). 

Let M, M' be any two matrices G SL(2,M), with parameter values 0,(j),X and 
0',(j)',X' respectively. Suppose that we make a retrodictively optimal, balanced 
measurement of the observables xm, Pm- Let fixM, Apa/ be the pointer observables 
representing the result of this measurement. Define new pointer observables /txA/', 

flPM'- 

/AxAfA ^/AxaA 

\^^PM' J \^J-PM J 

where L is the matrix 

ii-XM' J //pAf provide a measurement of x m' , Pm' ■ We now ask: what is the condition 
for this measurement to be retrodictively optimal and balanced, the same as the 
measurement of xm, PivH 

The retrodictive error operators for the M'-measurement are given by 

/ eXA/'iA ^ f ULcas AxAf'C/nicas ^ XM'\ ^ ^ ( ^XA/i 
\ePM'iJ WioasAPA/'C^mcas -Pa/'/ VPMi 

[see Eq. (|)]. Therefore 

i^XM'i) = Llxi^XMi) + I^lpi^PMi) + LxxLxp{{ixMi^PMi + epAfieXA/i)) 
i^XM'i) = ■^p2;(eXA/i) + ■^pp(epAfi) + LpxLpp({ixMiipMi + ipMi^XMi)) 

The fact that the M-measurement is retrodictively optimal and balanced means 
that (cxjv/i) — (^PA/i) — 5- Also, we have from Lemma 2, proved in ref. that 

(eXA/i + ^epAfi) IV" ® 0ap) = 

for every initial system state {ip) (where \(j>ap) is the initial apparatus state, as 
before). Hence 

{^(exA/iepAfi + epAfiCXA/i)^ = —i ^(exA/i + iepAfi)^^ = 
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Consequently 



i^XM'i) — 2 {-^xx + ^xp) 



i^XM'i) — 2 {-^px + Lpp) 

It follows that the Af'-measurement is retrodictively optimal and balanced if and 
only if 



which implies 



A- — A- =1 

XX ~^ xp px ^ pp 



■ LxxLpx ^' LxpLpp 1 



where L"^ denotes the transpose of L. Since detL — 1 we must have LxxLpx + 
LxpLpp = 0, which means that L is a rotation matrix. We conclude, that the neces- 
sary and sufficient condition for the M'-measurement to be retrodictively optimal 
and balanced is that 

M' = R^M 

for some [where is a rotation matrix as defined in Eq. (|^)]. If M and M' satisfy 
this condition we will say that the corresponding measurements are informationally 
equivalent, and we will write M ~ M'. 

If M ~ M' it means that a retrodictively optimal measurement of xg+^, Pe-<j> to 
accuracies ±^===5, ± . „, yields exactly the same information as a retrodic- 

V2sec20' Av'2soc20 

tively optimal measurement of xg/^^/, po'-^' to accuracies ^"^5=^1 ^2scc20' 
The condition for M M' can alternatively be written 

M'^^M' = M^M (11) 

In other words, two measurements are informationally equivalent if and only if the 
corresponding matrices define the same phase-space metric [see Eq. (j^)]. 

We conclude this section by showing, that given any M £ SL(2,R), it is possible 
to find Mq ^ M with zero obliquity (so that the axes are perpendicular). 

We can write 

M-M=(; 5) (12) 
where a, 6 > and ab — = 1. Let Mq be a matrix with zero obliquity: 

/ T- cos On sin 9q 
y— Ao sin 9q Xq cos 0oy 

[see Eq. (^)]. We want to show that it is possible to choose Aq, 6*0 so that 

M^Mo = M^M 

which means 



1 



^2 A^ ) cos 26*0 



a + b= { ^ + Xl 

1 
X 



2c = ( -2 - A;^ ) sin 29o 
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It is readily confirmed that these equations are soluble. An explicit solution is 

(13) 



^0 = 



itan-iy^ if 6^ a 

-4 if 6 = a 



and 

^((a + &)-sign(a-6)((a + &)2-4)^)^ if & ^ a ^^^^ 

\J a + c if = a 

Using these formulae, it is straightforward to express directly in terms of 

the parameters 6*, A. For the sake of simplicity we confine ourselves to the case 
A = 1, when 



70 = 




and 



Ao 



if 7^ 
if (7!) = 



(sec 20 + sign(sin 2ff) tan 2(/)) = if sin 26* 7^ 

(sec 20 + cos 26* tan 2(/))^ if sin 26*^0 



where the notation means "6* mod -I": 

^ 2 2 - 2 

for every integer n. 

This means that a retrodictively optimal measurement of the non-orthogonal 
quadratures ig+^j pe-^ to accuracies i"yf==5j and ^XTflocl^ yields the same 
information as a measurement of the orthogonal quadratures xq^ , pe^ to accuracies 
and zfc . Fig. ^ gives an illustration, for the case A = 1, = 0, </> = 40° 

[implying ^ -45° and Aq = (sec80° + tan 80°)^]. 

4. The Distribution of Measured Values 

Let M be any matrix S SL(2,R), and consider a balanced retrodictively opti- 
mal measurement of the observables XM^ Pm- Let flxM and flpM be the pointer 
observables describing the outcome of this measurement. Define fix, Ap by 

= M-i ( ™^ 
Ap / vApm / 

(c./. Eq. (|l0|)). Then Ax, Ap are the pointer observables for a measurement of x 
and p. This measurement will not be retrodictively optimal and balanced unless 
M'^M = 1. 

Let Pm be the probability density function describing the result of the measure- 
ment of xm, Pm] and let Qm be the probability density function describing the 
result of the measurement oi x, p. Then 

Qm{x,p) = Pm{x^i,Pj^j) 

for all X, p, where 



Pm 



M 
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Let and \{x,p)m) be the annihilation operator and (normalised) squeezed 
state defined by 

1 

3 



aj^j = — {xm + ipAi) 



aM\{x,p)j,.i) = {xM +ipM) \{x,p)j^,j) (15) 
Using the result proved in ref. j5| we have 

QaA^^p) = PaA^m^Pm) = -^{i^^P)M\p\i^^P)M) (16) 

where p is the density matrix describing the initial state of the system. We see from 
this that Qm is a generalised Husimi function of the kind defined in Section |l[ 
It was shown in the last section that there exists V' such that 



where Aq, are the quantities defined by Eqs. (13) and (|14|) , and where R^p, Tx^, 
are the matrices defined by Eq. (|^). In view of Eqs. (|6[|^) we then have 

\ix,p)^j) = e^^D,,Ulwlul\0) = e^>^a,,ulwl\0) (17) 

where x is a phase, D^p is the displacement operator 

D.j;p = cxp [-i{xM Pm - Pm Xm)] = cxp [-i{xp - pi)] 

and where |0) is the vacuum state annihilated hy a = -^{x + ip). Hence 

QMix,p) 

= ^ J dx'dp' exp [-a{x' - x)^ - 2c{x' - x){p' - p) - b{p' - p)^] W{x',p') 

(18) 

where a, b, c are the quantities defined in Eq. ( p^ and W is the Wigner function 
describing the initial system state. It can be seen that the quadratic form occurring 
in the argument of the exponential is the phase space metric corresponding to M 
[c.f. Eqs. (|) and (|l|)]. This means that two retrodictively optimal measurements 
define the same distribution if and only if they are informationally equivalent — as 
was to be expected. 

In terms of the parameters Aq, Oq we have 

It can be seen that Qm is obtained by smoothing the Wigner function on the scale 
Ao parallel to the X0„ axis, and on the scale parallel to the pg^ axis. 
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Figure 1. Oblique axes: parameterisation 




Figure 2. Equivalence of oblique and orthogonal phase space co- 
ordinate systems. Measuring the oblique coordinates 2:40 and p_4o 
each to the same accuracy i:yf=^gf = ±0.29 is equivalent to 
measuring the orthogonal coordinates X-45 and p-45 to accuracies 
= ±2.39 and ± = ±0.21 respectively. The measurement 
of p-45 is more accurate than the measurements of X40 and p-4o- 
The measurement of X-45 is much less accurate. 



